Generation of Entangled A/-Photon States in a Two-Mode Jaynes-Cummings Model 
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We describe a mathematical solution for the generation of entangled A'-photon states in two field modes. A 
simple and compact solution is presented for a two-mode Jaynes-Cummings model by combining the two field 
modes in a way that only one of the two resulting quasi-modes enters in the interaction term. The formalism 
developed is then applied to calculate various generation probabilities analytically. We show that entanglement, 
starting from an initial field and an atom in one defined state may be obtained in a single step. We also show that 
entanglement may be built up in the case of an empty cavity and excited atoms whose final states are detected, 
as well as in the case when the final states of the initially excited atoms are not detected. 
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I. INTRODUCTION 

Entangled states are one of the building blocks in quantum 
information processing and non-locality tests d. They can 
be used, in the case of the electromagnetic field, to improve 
the sensitivity of interferometric measurements HI Ell- 
and may help to overcome the classical Rayleigh diffraction 
limit in quantum optical lithography q|, A feasible way to 
generate such states is given by the atom-field interaction in 
the framework of one- or two-mode Jaynes-Cummings (JC) 
models 

Here we consider the generation of entangled two-mode 
field states by different schemes inspired partly by Refs. [7;, 
l8l fl3ll . We let two-level atoms interact, one at a time, with 
two degenerate modes of a lossless cavity. Solving the corre- 
sponding JC model algebraically by an SU(2) transformation, 
we discuss the generation of entangled -photon states of the 
general form 
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which comprises the maximally entangled Bell states 



|^> = ^=(|JV,0)±|0,AO). 



(2) 



II. ALGEBRAIC SOLUTION OF THE TWO-MODE 
JAYNES-CUMMINGS MODEL 



The JC Hamiltonian for resonant interaction of a two-level 
atom (\e), \g)) with two field modes (a\,a 2 ) in the dipole and 
rotating wave approximation is given by H = Ho + H ml , where 
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Herec,:^ \e)(e\-\g)(g\,G+ := \e)(g\,G := \g)(e\ and 1 = 
\e)(e\ + \g)(g\ are operators for the two-level atom, gj is the 
coupling constant of the ;th mode with the atom, and fm is 
the photon energy. We introduce the quasi-mode operators 



M = Yifli +Y2fl2, A 2 = -f 2 a\ + Yi fl 2, 



(5) 



where y ; := g t /g, and g := VL?i I 2 + 



Equation (0 de- 
fines an SU(2) transformation of the mode operators a\,a.2, 
leaving the commutation relations and the number-sum oper- 
ator a[a\ +a\a 2 = A\A\ +Ai,A2 invariant. The transformed 



Hamiltonian then reads 



The field states are defined in terms of the usual two-mode 
Fock states \n\,n 2 ) := \n\)\\n 2 ) 2 , with n\ {n 2 ) photons in 
mode one (two). The two modes have the same energy and are 
in resonance with the two-level atom. We solve the model al- 
gebraically by combining the two field modes into two quasi- 
modes of which only one enters in the interaction term, yield- 
ing an effective one-mode JC model 0. Using its known so- 
lution and the transformation between mode and quasi-mode 
Fock states, the generation probabilities of the entangled states 
are found for three different schemes. 
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representing a JC Hamiltonian for the quasi-mode A i decou- 
pled from the non-interacting quasi-mode A 2 . Since H^ t de- 
pends only on quasi-mode one and [//o,//i n t] = 0, the time 
evolution operator U(t) = exp(— iHi nt t/h) in the interaction 
picture is the same as for a one-mode JC model. Expanding 
U in the atom basis {|e), 
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U = U ee \e) (e\ + U ge \g) (e\ + U eg \e) (g\ + U gg \g) (g\, (8) 
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the matrix elements U a b{t) are given by [15| 
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where T := gt is the dimensionless "interaction time". The 
model can be solved in the usual way in terms of quasi-mode 
Fock states defined as the common eigenstates of AjAi and 
A\A 2 . The complete solution is then found by giving the rela- 
tion between the quasi-mode and the mode Fock states. 

The quasi-mode operators A,-, A], i = 1,2, obey the same 
algebra as the mode operators a,-, a\, so that two-quasi-mode 
Fock states (denoted by a double-ket) can be defined by 
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To find the transformation between the two-mode Fock states 
\n\,n 2 ) and the two-quasi-mode Fock states \n\,n 2 )), we use 
Schwinger's oscillator model Jl6ll and introduce angular mo- 
mentum states \j,m) and \j,m)), where j = {n\ +n2)/2 and 
m = (n\ — n 2 )/2. In cases where it is not obvious, we 
shall write a subindex S on the state vectors to indicate the 
Schwinger angular momentum basis, e.g., |2,0) = 1 1 , 1 ) s - In- 
serting Eq. (|5} into Eq. ilOi and identifying the two vacua 
|0,0» and |0,0), we obtain 



\J,m)) ■ = 



f x a\+f 2 a\ 
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Expanding the products, rearranging the terms ifTvtl and using 
the definition of the Fock basis \n\ ,n 2 ) in terms of a\ and a 2 , 
we obtain the important relation between the quasi-mode and 
the mode Fock bases 



\j,m)) = £ D$ m to,*,x)\jM), (ID 

m'=-j 

\j,m) = £ D { £(^, X )\j,m')). (12) 



The action of U a b on the field states is easily calculated in 
the quasi-mode Fock basis 

U ee {%) \j,m}) = cos (xy/j + m+l) |j, m )), 

U ge (x)\j,m}) = -isin (xyj j + m + l) \ j + \.m + ±)) , 

U eg (x) \j,m)) = -isin (xy/j + mj \j -\,m- ±» , 

U gg (T:)\j,m)) = cos (xy/J+mj \j,m)), (13) 



showing that U ee and U gg do not change the number of quasi- 
photons, whereas U ge (U eg ) act as creation (annihilation) op- 
erators of quasi-mode one. Using Eq. il It and Eq. d!2i . we 
find for the action on the usual Fock states 
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where we have introduced the following coefficients 
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Given the above equations, we now have all the ingredients to 
calculate the time evolution of the density operator according 
top(t) = U(t)p(0)U f (t). 



III. GENERATION OF ENTANGLEMENT IN ONE STEP 



Here D if»(«P. lC |.X) = exp[-i(m'9 + m%))d^Jp) are 4e 
Wigner D-matrix elements of the SU(2) group lUfl fl7ll . 
with arguments determined by (p = (pi — 92, X = <Pi + <?2, 
cos(^/2) := |yi|, sin(#/2) := |y 2 |, and y = |y,|exp(icp,). It 
follows that the mode and quasi-mode Fock states belonging 
to the same total number of photons, n\ + n 2 = 2 j, are related 
by an irreducible rotation matrix of weight j and with Euler 
angles determined solely by the interaction constants. 



We start with the calculation of the probability to find at 
time X the field state |\PV) in Eq. Q, assuming an initial field 
state |^) and an atom entering the cavity in either the excited 
or ground state. The analytical calculation is straightforward. 
The initial field state is expanded according to 

= E E b nm 2 \ni,n 2 ) = £ £ fym !./',«), (16) 

n 1 — ri2=0 j=0 m= — j 
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where the primed summation symbol indicates a sum over in- 
teger and half integer values of j. The expansion coefficients 
with respect to the Fock and Schwinger basis are related by 



u J+m,j-m 

basis by 



= b 



j.ni- 



The state I^n) is given in the Schwinger 
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cn 
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with density operator p^ N = ^^(IVl. From the time- 
evolved initial states 

U\e;%) = U ee \e-£)+U ge \g-£), 
U\g-£) = U es \e-£)+U gg \g-£), 

we obtain the reduced density operator of the field by 
tracing out the atomic degrees of freedom: pp (0 = 
tr A (U(t)\a;t,)(a;^\U f (t)), a = e org. The probability to find 

I'Pjv) at time t follows from (p$) = tr(p^\t)p\y N ) and is 
given by 
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for the initial atom- field state and by 
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(19) 



for the initial state |g; ^) . It follows that in order to obtain non- 
vanishing probabilities at time x, the initial field state must 
contain at least one of the Fock states from the set 

{|Af,0),|jV-l,l),...,|0,AO} 

U{|tf-l,0>,|tf-2,l>,...,|0,tf-l>}, (20) 

if the atom is initially in the excited state, or from the set 

{|JV,0),|jV-l,l),...,|0,AO} 

U{|tf+l,0>,|tf,l>,...,|0,tf+l>} ) (21) 

if it is in the ground state. 

For = 1 the set of contributing atom-field states according 
to Eq. d20b is 

{|e;l,0),|e;0,l)}U{|e;0,0)} (22) 
and according to Eq. i !2 It 

{|g;l,0),|g;0,l)}U{|g;2,0),|g;l,l},|g;0,2)}. (23) 




FIG. 1: Parametric plot of the generation probabilities (p<j/+) (solid) 
and (pip-) (dashed) as function of time x = gt, for different initial 
atom-field states shown at the right. 



We illustrate the probabilities for the generation of l^f) = 
(|1,0) ± |0, 1))/V2 in Fig.^ where we have taken gi = g2, 
91=92 = (real coupling constants). The states shown on 
the right are just the initial atom-field states from Eqs. (122b 
and J23I . The interesting case is |e;0,0), where the state 
is produced periodically with probability one at the times X„ = 
(n + l/2)/jt,forn = 0,l,... 

Next we consider the creation of various Bell states, 
Eq. 10, resulting from the initial atom-field state \e;N,Q) = 
l e ' T' t)s - ^ e orJta i n fr° m Eq- < l 1 t>l > the entangled Af -photon 



field states j 1 ?^) with probabilities 



Pm± 



cj N (x)±ci 
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,(x) 



(24) 



as well as the entangled (N - 
probabilities 



1) -photon states with 



T W4 



N 

° N±l N 
T 



(T)±S' 



,(x) 



(25) 



In the case of Eqs. d25i the Bell states l^^j ) have no overlap 
with the initial field state \N, 0) . The probabilities at time x, 
however, may come close to one for some particular values of 
the coupling constants and interaction time. In this case we 
may say that |M*jy +1 ) have been generated in a single step or 
single shot. This property of the two-mode JC model can be 
understood if we think of the atom (re)emitting photons into 
and (re)absorbing photons from the two modes many times 
during the interaction time x. 



IV. CONDITIONAL GENERATION 

Next we present a conditional scheme for the generation 
of jV-photon entangled states starting with an empty cavity 
010] ■ The scheme implies sending consecutively atoms in the 
excited state through a two-mode cavity and detecting them 
in the ground state. We start with an initial atom-field state 
|e;0,0) = |e;0,0)) s and let the first atom interact for a time 



Xi. By using Eq. d 1 3I > we obtain the state 



t/(xi) k;0,( 



; cos(xi) |e;0,( 



isin (xi) 



2' 2))s ' 
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Detecting the atom in the ground state leaves the field 
in the state = K\(— i)sin(xi) |j,j)) s . where K\ = 
| sin (xi ) | ~ 1 exp(iai ) is a normalization constant. By choosing 
the phase 0Ci appropriately, the factor entering the normalized 
state may be set equal to one, yielding the state 1 5 , 3 )) s " Pro- 
ceeding this way the field state obtained after N conditional 
steps is simply given by 



\%n) 



IN N 
I 2 ' 2 



))s=L D N: k N(<?,®,X)\N-k,k), (26) 



k=0 



where we have used Eq. il It and Fock-state notation on the 
r.h.s. This is precisely a state of the form given in Eq. Q 
with coefficients determined by the Wigner rotation matrix el- 
ements. Since these elements depend solely on the coupling 
constants, the generated entangled state is sensitive to their 
magnitudes and phases. The state in Eq. i2bt corresponds 
to the quasi-mode Fock state \N,0)), implying that each con- 
ditional step generates one photon in quasi-mode one. The 
generation probabilities of the states 
conditional steps are given by 



and |^) after N 



\(Vn\Xn)\ 2 
\(^n\Xn)\ 2 



N/2 
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N N zszu N N 
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(27) 
(28) 



We shall show that the probability to detect the atoms N times 
consecutively in the ground state is a rapidly decaying func- 
tion of N. But, as discussed below, it is not essential to rely 
on this assumption. Actually, it is sufficient to detect them in 
a sequence of n(> N) steps N times in the ground state. 



V. NON-CONDITIONAL GENERATION 

In the following we consider a non- conditional scheme. We 
start with an empty cavity and send a sequence of excited 
atoms through it without detecting their final states. The re- 
duced density operator of the field after the passage of the first 
atom (interaction time Xi) is given by 



p f 1} (xi; 



cos 2 (xi) |0,0)> s ((0,0|+sin 2 (xi) \\,k)) s ih\ 



and serves as the "initial" field configuration for the second 
excited atom. Proceeding this way, the reduced density oper- 
ator of the field after « steps turns out to be of the form 

n/2 

PF ({*«})=£' pf({*n})\j,j)) S (UJl (29) 

;=o 



where the coefficients p" are given recursively by 

M 2i \ (n-l) 

p [ Q > = cos z (x„)^ ; , 
pf = cos 2 (x n V27TT)pf- 1) + sin 2 (x„ 



2j) PjA , 



_(») 

Pn/2 



sin 



(«-l) 
l)/2' 



(30) 



for 1/2 < j < (n — l)/2 and p^ = 1. The argument {x n } 
stands for all interaction times (xi,...,x„) of the n steps. 
Equation ( I29> . which is obviously true for n = 1 and n = 2 
(see Eq. OH ), can be proven by induction. 

The coefficients pf in Eq. i29i are the probabilities to 
find the field after « non-conditional steps in the quasi-mode 

state |y',j))s- In particular p^ = cos 2 (xi)cos 2 (Z2V2) . . . 
cos 2 (x„ v / n) and pf 2 = sin 2 (xi)sin 2 (X2v / 2) • ■ -sin 2 (%„y/n) 
correspond to the cases where in n steps the initially excited 
atoms emerge n times in the excited and ground state, respec- 
tively. The intermediate pf's correspond to the cases where 
the n atoms emerge 2 j times in the ground state and n — 2j 
times in the excited state, irrespective of the order of appear- 
ance. The coefficient pf consists of a sum of (£.) terms, 
each of which corresponds to a particular sequence of |g)'s 
and |e)'s contributing, respectively, a sine squared and cosine 

squared factor. There are altogether 2" terms in Eq. ( I29l l. All 

(2) 

this is easily seen by giving Pp as an example: 



p(f) = C0S 2 Xi COS 2 X2 |0 



j a g {(0, 

2^ „;„2„ , „;„2„ „„„2 



cos Xi sin 12 + sin i\ cos (X2V2 
+ sin 2 Xisin 2 (x 2 V2) |1, l)) s ((l, 1 1 . 



I 2 ' 2))s(( 2 1 2 I 

(3D 



Here the four terms correspond to the final state sequences 
(e,e), (e,g), (g,e), and (g,g). 

The states 1^) and |^) are generated in a non- 
conditional 77-step process with probabilities 



P>i'= 



P N 
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(32) 
(33) 



which are the conditional probabilities found before, multi- 

(n) 

plied by the probability p N j 2 - H ere lh e interaction times must 

In) 

be chosen such that p N j 2 0, which amounts to control the n 
parameters (xi , . . . ,x„). The state ^n) can be generated in a 
minimum number of n = N steps with probability p^ 2 which, 
however is a rapidly decaying function of Af. 

In the non-conditional scheme all field states \j,j)) s = 
\2j, 0}) for j = 0,l/2,...,rt/2 are produced, Eq. <|55}. On 
the contrary, in the conditional scheme only the entangled Af- 
photon state |Af,0)), Eq. &26\ . is generated, if in n steps Af 
atoms are detected in the ground state. To produce it 
is, therefore, not crucial that the atoms have been detected Af 
times consecutively in their ground state. Any sequence of 
ground and excited states containing Af times the ground state 
will do it. Finally, we note that there is a particular choice of 
the interaction time of the £th atom, given by z e = n/(2V£) 
for which both, the conditional and non-conditional scheme 
give (with probability one) the same entangled state |Af , 0)) in 
Eq. (El. 
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VI. CONCLUSION 

To conclude, we have solved the two-mode JC model al- 
gebraically by reducing it to an effective one quasi-mode JC 
model. The mode and quasi-mode picture are unitarily re- 
lated by an SU(2) transformation. The solution found is used 
to discuss three different schemes for the generation of en- 
tangled states of the two field modes. To generate entangled 

-photon states in a single step the initial state must contain 
at least N — I photons and an excited atom. Starting from the 



vacuum we need at least n>N steps to produce pure (mixed) 
field states in the conditional (non-conditional) scheme pre- 
sented. 
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